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The relativistic theory of above-threshold ionization (ATI) of hydrogen-like atoms in ultrastrong 
radiation fields, taking into account the photoelectron induced rescattering in the continuum spec- 
trum is developed. It is shown that the contribution of the latter in the multiphoton ionization 
probability even in the Born approximation by Coulomb field is of the order of ATI probability in 
the scope of Keldysh-Faisal-Reiss ansatz. 



I. INTRODUCTION 



The increasing interest to the process of multiphoton above threshold ionization (ATI) of atoms in superintense 
laser fields particularly is conditioned by the problem of high harmonic generation and short wave coherent radiation 
^ implementation via multiphoton bound-free transitions through free continuum spectrum (one of the possible version 
of a free electron X-ray laser). During the last two decades numerous investigations have been carried out to study 
ATI of atoms both theoretically and experimentally and many review papers (see, e.g., |l0|| ) and monographs 
pT|- are devoted to this problem. 
The main ansatz in the nonrelativistic theory of multiphoton ionization of atoms in strong electromagnetic (EM) 



' radiation fields is the Keldysh one |16| (further called Keldysh-Faisal-Reiss ansatz IT^, [H). The advantage of this 



0^ [ approach is that it leads in a very simple way to reveal some of the main qualitative features of the photoelectron energy 
. spectrum in ATI experiments ( jl^ , and ||2l[ ) . Within the scope of this ansatz the photoelectron rescattering in the 
O I field of atomic remainder is neglected. In order to cover this gap attempts have been made to describe the photoelectron 
^ ' , , final state by " Coulomb- Volkov" wave function that is a product of the Coulomb wave function of clastic scattering 
^ ■ and a wave function of electron in the EM wave field This wave function results in the factorization of the 

Ci probability of multiphoton ionization and restricts both the frequency (low frequency approximation) and intensity of 
^ , the wave. The use of another ansatz for definition of multiphoton ionization probabilities ps| ] should also be noted. 

The description of the photoelectron final state taking into account the stimulated bremsstrahlung (SB) at the 
photoelectron scattering on the electrostatic potential of the ionized atom in the presence of strong EM radiation field 
(induced free-free transitions) still remains as one of the main problems for the ATI process. Moreover, the definition 
of dynamic wave function of an electron in SB process already is problematic, so that the main results concerning to 
multiphoton SB probabilities have been found trough S-matrix formalism for " free- free" (over electrostatic potential) 
transitions in the Born approximation between Volkov states in EM wave pO[ |. Although, in a lot of cases when 
the condition of the Born approximation is broken, the scattering process is described in low-frequency p2[ or 
eikonal ||3^ approximations. Though the Born and low-frequency approximations in SB process are applicable for 
describing free-free transitions in high intensity radiation fields, but they don't take into account the mutual influence 
of the scattering and the radiation fields [i.e. the probability of SB is factorized by elastic scattering and photon 
emission or absorption processes]. What concerns to the eikonal approximation in SB process it is not applicable 
beyond the interaction region [ z ^ pa? /Ti, where z is the coordinate along the direction of initial momentum ~p* of 
the particle, a is the range of the interaction region, and Ti is the Plank constant]. The description of the electron 
eigenstates in SB process behind of the scope of these approximations has been made in , developing a generalized 
eikonal approximation (GEA). The obtained GEA wave function enables to leave the framework of ordinary eikonal 
approximation and abandons from the restriction z <C pa? /%. Besides, such wave function simultaneously takes 
into account the influence of both the scattering and radiation fields on the particle state. So, to determine the 
multiphoton probabilities of above threshold ionization of an atom it should be known the wave function of the 
ejected photoelectron in SB process with more accuracy. On the other hand, in the current superintense laser fields 
the state of electron becomes relativistic already at the distances Z <C A (A is the wavelength of a laser radiation) 
independent on its initial state. Hence, the problem of ATI of atoms with the photodetached electron SB process 
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demands a relativistic consideration. 

The relativistic generalization of multiphoton SB in the first Born and eikonal approximations have been made 
in the papers and respectively. In |Q on the base of the solution of the Dirac equation the GEA 

approximation ]34| has been developed for relativistic scattering theory in the arbitrary electrostatic and plane EM 
wave fields, including both the Born and eikonal approximations in corresponding limits and describing the spin 
interaction as well. Such a wave function allows us to describe the final state of the photoelectron with more accuracy 
in the ATI process of atoms. 

The relativistic description of multiphoton ATI of hydrogen-like atoms for high-intensity laser fields taking into 
account the spin interaction has been developed analytically in the papers |3^- with an approximation where the 
stimulated bremsstrahlung of the emergent electron is neglected. The relativistic consideration of ATI is important 
as it is generally assumed that the problem of stabilization of atoms in ultraintense laser fields must be solved within 
the framework of relativistic theory 1 4^ . From this point of view some attempts have been made to solve analytically 
the Klein-Gordon equation or numerically the Dirac equation Q|, in fields of a static potential and 

monochromatic EM wave [using various model potentials of one or two dimensions and various approximations]. In 
the papers [Q- the relativistic corrections to the nonrelativistic results have been given. 

Note, that at the present time an analytic formulas for these probabilities, taking into account the photoelectron 
rescattering, are unknown even in the first Born approximation for the Coulomb scattering field. So, in the present 
paper the relativistic probabilities of multiphoton ATI in the limit of the Born approximation for the photoelectron 
rescattering, are calculated. Moreover, it is shown that the neglect of the photoelectron rescattering in the relativistic 
domain specially |^^- is invalid, since the contribution of the electron rescattering process in the matrix elements 
of transitions has the same order by a scattering potential in the Born approximation as the matrix elements of 
bound-free transitions for the ATI process. 

The organization of the paper is as follows. In Sec. II we present the multiphoton cross sections of the above 
threshold ionization of hydrogen- like atom in ultraintense laser field (with the help of GEA wave function), taking 
into account the induced free- free transitions of the ejected photoelectron in the continuum spectrum. Because of much 
complicated expressions in GEA the spin interaction is neglected and the ultimate analytic results for the multiphoton 
probabilities are performed in the limit of the first Born approximation by ion (atomic remainder) potential which 
we present in Sec. III. In Sec. IV we treat the dependence of ATI probability on polarization of an EM wave and 
consider the differencies between circular and linear polarizations of electromagnetic wave. 



II. THE IONIZATION PROBABILITY BY RELATIVISTIC GEA SOLUTION OF A WAVE EQUATION 

OF AN ELECTRON 

The problem has been reduced to the investigation of the relativistic exploration of the transition S-matrix formalism 
utilizing the relativistic GEA wave function ||38|| as a wave function of the final state of a photodetached electron (it 
has been neglected with the spin interaction in relativistic GEA wave function). Following the relativistic S-matrix 
formalism the bound- free transition amplitude can be written in this integral form (in natural units h = c = 1) 

/oo 
^'(-)t (x)V<i>{x)d'^x , (2.1) 
-oo 

where x — (i, ~r*) is the four-component radius- vector x^, ^{x) is the initial unperturbed bound state of the atomic 
system and 5''^~^(x) is the final out-state of an electron in the potential of atomic remainder and in the field of a plane 
EM wave {K'' is denotes the complex conjugation of K). We assume the EM wave to be quasimonochromatic and of 
an arbitrary polarization with the vector potential 

A ((ys) = Ao(iy9) (l?! cos -|- if 2C sin (/?) ; ip — k ■ x ^ cut — fc -T*, (2-2) 

where k — (oj, fc ) is the four- wave vector, ^0(1^9) is the slow varying amplitude of the vector potential of a plane wave, 
"e*i and "6*2 are unit vectors: ifi _L 1?2 -L k , and arctan^ is the polarization angle. 
According to the Klein-Gordon equation the interaction operator is 

V ^-1e~l{<f){-i^)^e^~t{Lp), (2.3) 

where e is the electron charge. 

The wave function of the final state of the photodetached electron in the relativistic GEA approximation has the 
following form ( [|8|) 
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{x) = -^F\x)exp[-tSv{x)], 

V ^i^O 



The Svix) is the action of photoelectron in the field (2.2) 



Sv{x) = Tf-T* -Hoi + a 



^ ^ shi[ip - 0it)] - ^(1 - C) sin 2^. 



(2.4) 



(2.5) 



Here 11 — (XIq, T? ) is the average four- kinetic momentum or " quasimomentum" of the electron in the plane EM wave 
field, which is defining via free electron four- momentum p — (eg, ~jf) and relative parameter of the wave intensity Z 
by the following equation 



n = p + kz{i + c^y, 



z = 



(2.6) 



where is the averaged value of the amplitude Ao{ip). The wave function (2.4) is normalized for the one particle in 
the unit volume V — I. 

Including in (2^) quantity (^j^^ is the intensity-dependent amplitude of the electron- wave interaction and as a 
function on any three-vector b has the following definition 



with the phase angle 



~jt) = arctan I ^ C, 



(2.7) 



(2.8) 



The function F^{x) in Eq. (^!4|), which describes the impact of both the scattering and EM radiation fields on the 
photoelectron state simultaneously, has the following form 138] 



F^{x) 



exp 



4-^-3 



E 



[a (g) Dl^ {0,{lt) 9ijt)) Z{1 e)Dl^] UpDj} 



+ 2Tf- Y - 2n(fc • p - + «0 



xU{~q') exp [—i {"^ ■ + ai ( sin [(p — 9i ( Y)] — "2 ("g^) sin 2ip + 6i {~^)n} 



where 



U{lf) = / [/(T*) exp(— «"g*-T^)(iT* 



(2.9) 



(2.10) 



is the Fourier transform of the potential of the atomic remainder and ai{lf) , a2{~<f) are dynamic parameters of the 
interaction defining by expression 



and Oi{lf) is the phase angle 



k ■ if 



2{k ■ p ~ k - if 



(2.11) 



e,{if)^e(^it ■if)r/k-p + ify (2.12) 

The functions J„(w, v, A), _D„, Di ^ (^i(~9*) — , and 132, n are defined by the expressions (also see Ref. ]3^ ]) 

Dn - J„(ai( q^),-a2(Y),0i( q^)), (2.13) 



D^,^{e^{-7i)-e{-f))^- 



J„_i(ai(^),-a2(Y),ei(^))e-<«^(^)-^(f^') 
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+J„+i(ai(^), -asm, 0i(Y))e<'^(^'-'(^)) 



(2.14) 



and 



Do 



Jn.2ia,ilf), -a2ilf), 0i(Y))e-^2ei(^) 



(2.15) 



In the denominator of the integral in expression (2.9) +iO is an imaginary infinitesimal, which shows how the path 
around the pole in the integrand should be chosen to obtain a certain asymptotic behavior of the wave function, 
i.e., the ingoing spherical wave [to determine that one must be passed to the limit of the Born approximation at 

Since we consider the ATI problem for hydrogen-like atoms ( Za ^ 137), the initia l ve locities of atomic electrons 
are nonrelativistic and as a initial-state wave function <f> in the transition amplitude (|]^) will be taken a stationary 
wave function of hydrogen-like atom bound state in nonrelativistic limit 



$(T*,t) 



/2m 



$0 (^) cxp {-isat) 



Eq = m — El 



where Eb )^ is the binding energy of the valence electron in the atom 

2mEB = a^^. 



(2.16) 



(2.17) 



Concerning the relativism of the photoelectron final state in strong EM field it is followed to mention that at the 
wave intensities already ^ 10""'^, where 



eAp 
m 



(2.18) 



is the relativistic invariant parameter of the wave intensity, relativistic effects become observable and the final state 
of the photoelectron should be described in the scope of relativistic theory. Moreover, at the currently available laser 
intensities C ;^ 1 (even ^ ^ 1) a free electron becomes essentially relativistic already at the distances smaller than 
one wavelength. On the other hand, in such fields becomes actual the production of electron-positron pairs from 
intense photon field on the electrostatic potential of atomic remainder through multiphoton channels. However, we 
can calculate separately the ATI probability in superstrong laser fields without restricting intensities by the threshold 
value of multiphoton pairs production (^ ~ 2; see and ^|) since those are independent processes. 
Since y is a Hermitian operator, the transition amplitude (2.1|) can be written in the form 



T,. 



<i>{x)V'' {x)¥-^ ''{x)d'^x. 



To integrate this expression it is convenient to turn from variables t, to ip, if [see (2.2)] 



T, 



and make a Fourier transformation of the function (x) over variable ip 

oo 

FHv'^V')^ Y1 ^;(^)exp(-iZ(/3), 



(2.19) 



(2.20) 



(2.21) 



/— — OO 



1 

2^ 



F{ip, if) ejip(ilLp)d(p. 



(2.22) 



Then with the help of Eqs. (p|)-(pl|), ( pl6| ) [using as well Eq. ([A^)] and taking into account the Lorentz condition 
for the plane wave field k ■ A{ip) — 0, we can acco mpli sh the integration over the variable f in Eq. ( 2.2C| ). After a 
simple transformation with the help of the formula (A5) we obtain the next expression for the transition amplitude 
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L.I — — O0 



k ■ p 



xd L - I 



+2 E y (4^^'(^+^)^(^ 



xa 



[lu [a (g) i?i_„ ie.it) " 0( P^)) - Z(l - e)Dln] - UoDi} 



2rf - 2n{k ■ p - ~t-~ct) + iO 



, Ho — eo r ; , 

X L — I + n 



(2.23) 



where ~^ is the three- vector 



9 = P 



(2.24) 



and the function <i>; ^ 6 ^ is the Fourier transform of $i(lf) = $(lf )F;("^) and as a function on any three- vector b is 
defined by the formula ( 2.10 ), and 



J„„i(ai(7f + 9^), -f (1 - C%OiTf + Y))e-<^(P^+^)-«(^)) 



-J„+i(oi( g^), -f (1 - C^), 0{-jt + if))e<'(T^+^)-^(^)) 



(2.25) 



where the parameters a ( — Ejt^ — \ 0(jjf + if) are determined by the expressions (2.7) and (2 
\fe p- k ■~q^ ' — 



and 



—r2 



4(fc ■ p — k ■ if) 



(2.26) 



Using the general conservation law of considering process the probability amplitude of the above-threshold ionization 
in concluding form can be presented in this ultimate form 



T, 



i27r(fc • p) 
\frn 



'w^ E {(A^ + ^-^(i + c'))$K: 



N,l=- 



5 



xa 



+i(Ar_i+„)e(p'+"g') 



"g*^ + 2rf-"g* - 2n{k -p- k -it) + iO 



X 5 (IIo - eo - LoN) . 



(2.27) 



The differential probability of ATI process per unit time in the phase space dlt / {2Tr)^ (space volume F = 1 in 
accordance with normalization of electron wave function) taking into account the all final states of photoelectron with 
quasimomenta in the interval Tf , Tf + dlt is 

dli 

dWt^f = Wi^f- — -o- 
{27rf 



dn 

(2^' 



where dft is the differential solid angle and 



is the "effective mass" of the relativistic electron in the EM wave field. 



2{i + e) 



Including in Eg. ( 2.28 ) the transition probability per unit time Wi^f has the general definition 



Wi^f = lim i \T,^f\^ 

i— too t 



(2.28) 



(2.29) 



(2.30) 



III. THE RELATIVISTIC BORN APPROXIMATION BY THE POTENTIAL OF ATOMIC REMAINDER 

FOR HYDROGEN-LIKE ATOM IONIZATION 

The impact of rescattering effect on the ATI process is more transparent in the limit of the Born approximation 
by the scattering potential. The latter takes place if the corresponding part of the action 5*1 (~f^,i) in the GEA wave 
function (2.4), i.e., the exponent of the function F^'lx) (F^ (x) = exp [iSi{'r^ ,t)]) is enough small (see Ref. [p8[): 



1^1 i^,t)\ « 1 



(3.1) 



Expanding Eq. (2.27) into the se ries and keep ing only the terms to the first order over J7(T*), after a simple 
transformation, utilizing Eqs. (A2), ( [A^ ) and (|A5|), we obtain 

■9 °° r ^ 



"v/mllo 



N=-oo 



.,^),-f(l-a^(l?) 



n— — 00 



(2^ 



[iv + ,^-Zi(i + c')] (k-p-1^ ■it]^{-^ + ii)u{it) 



6 



xe 



-in0i( q )+i{N+n)e( p + q) 



uj < a 



£>t,„(9i(t) - 9(ff)) - z(i - c')dL.} - noD;, 



■q'^ + 2rf --^ - 2n{k ■ p - + iO 



X 5 (Ho - eo - i^N) . 



(3.2) 



For hydrogen-like atoms with the charge number Za the condition of the Born approximation for the photoelectron 
scattering (in Coufomb field) 



hv 



< 1 



(3.3) 



requires electron velocities v » ZaU, where a = /Tic — 1/137 is the fine-structure constant (it is assumed that 
Za « 137, fi and c are restored for clarity). The photoelectron acquires such velocities in the EM wave field at the 
intensities 



Za_ 

137' 



(3.4) 



As will be shown bellow the Eq. ( |3.4D is the condition of the Born approximation in ATI process of hydrogen-like 
atoms taking into account the photoelectron rescattering. 



The initial bound state enters into Eq. (2.27) through its momentum space wave function $ (^ j . For hydrogen-like 
atom the bound state wave function is 



exp(-77/a) 



(3.5) 



where a — a^jZa (ao is the Bohr radius) (3.1) and the corresponding momentum space wave function has 

the following form: 



23(7ra3)i/2 



— >4 , 

b a4 



(3.6) 



Note, that in Eq. ( |3.6| ) it has been taken into account that b a 3> 1 in accordance with the Born approximation. 

Then the function $ + if) in the second term in figured brackets can be replaced by the quantity 5{lf + ~(f) /Vtw? 
because of the small contributions of the other terms in expansion of Ti^f over parameter ~^'^a? [see, e.g., [5l| ]] which 
will be shown bellow. Such a delta function can be used to accomplish the integration over if in the second term of 
the sum in the figured brackets of Eq. ( |3.2D . 

For the scattering of a charged particle in the Coulomb field for which the Fourier transform is 



An 



am g ' 



(3.7) 



we have the following expression for the transition amplitude in the field of arbitrary polarization of EM wave 

p \ Z. 



N=-c 



k-pj ' 2 



^£o 9 ' 
m{k-p) J: 



£ (2n-a'(l + C^)) 
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ml+sl~2eo{Tlo+Lo{2n~ N)) 



S (Uo- So - LoN) , 



(3.8) 



where a' is defined by Eq. ( 2.26| ) at ~ — "g^and a' = e'^AQ/Aujeo , then J„ " •' ^ is the ordinary Bessel 
function [ J2„ (0, a;,0) — Jn{x) (A7)], Cg and C2,s are defined by the expressions 



= Js{a{Tt/kp)), {Z - a'){l - e)/2, 9{lf)), 



and 



1 r 



2,s 



Js-2{a(jf/kp)), [Z - a'){l ~ C')/2, 0(l?)e-"''( P ) 



(3.9) 



+ J.+2(a(l?/M), - «')(! - C')/2, 0(l?)e^2^( P ^ 



(3.10) 



To determine the differential probability of ATI it should be integrate the expression ( |3.8| ) over Ho according to 
( ^.28 ). In result we have 



dW, 



dfl 



24 - (iv-z(i + c2))'(fc.n)2 



N=No 



2m{N - Z{l + C)){k-U) ^ 



where ^ is the three- vector 



((£0 + 2na;)cL2„ + - C')C^2V-2„) } 



= If - N^, rf = v/(eo+t^A^)2-TO2 



(3.11) 



(3.12) 



The number Nq over which is carried out summation in Eq. ( ^.ll| ) is defined from the energy conservation law of ATI 
process: A^o = {{m^ — eo)/uj) 



The first term in the figured brackets of Eq. ( 3.11 ) corresponds to result of the KFR approximation. And the 
second term shows the dependance of ATI probability on the ejected photoelectron stimulated bremsstrahlung (SB) 
probability, i.e., it takes into account the rescattering process. 



IV. PROBABILITY OF ATI PROCESS FOR THE CIRCULAR AND LINEAR POLARIZATION OF EM 

WAVE 



The state of photoelectron in the field of a strong EM wave and consequently the ionization probability essentially 
depends on the polarization of the wave [nonlinear effect of intensity conditioned by the impact of strong magnetic 
field]. Thus, for the circular polarization the relativistic parameter of the wave intensity = const = and the 
longitudinal velocity of the electron in the wave vn = const (eliminating this inertial motion - in the framework 
connecting with the electron - we have the uniform rotation in the polarization plane with the wave frequency 
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w), meanwhile for the Hnear one = cos^ </? and vn osciUates with the frequencies of all wave harmonics nu 
corresponding to strongly unharmonic oscillatory motion of photoelectron. The later leads to principally different 
behavior of ionization process and corresponding formulas depending on the polarization of strong wave. Therefore, 
we shall consider t he c ases of circular and linear polarizations of EM wave field separately. 

From the Eq. ( [j.ll| ) for the circularly polarized wave = 1) in the first Born approximation by ionized atom 
potential we obtain the next formula for differential probability of ATI process 



E 

N=No 



{N~2Zf [k ■ n) 



if 



72 



a 



1 



9 ■ 



2 (A^ - 2Z) {k ■ n) 



(4.1) 



As is seen from this formula in contrast to the case of another polarizations the differential probability of ATI 
process is defined by ordinary Bessel function instead of the function J„(m, u,A) and the sum over n vanishes. 
The latter corresponds the above mentioned fact that for the circular polarization the parameter of intensity of the 
wave = const and effect of intensity of strong wave is appeared in the form of constant renormalization of the 
characteristic parameters of the interacting system. 

Let us estimate the contribution of photoelectron rescattering in the probability of ATI process that is the second 
term in the figured brackets in Eq. (4.1). The latter is 



2{N-2Z) [k ■ n) 



1 



(4.2) 



for the most probable number of absorbed photons at which the Bessel function has the maximum value. So, 
the rescattering effect has the same order as the probability of the direct transition in SFA. It is followed to note 
that the derivations rely upon the SFA (e.g., |^1|]) are expected to become more accurate at high intensity EM 
field. However, the prediction of SFA regarding to rescattering effect in high intensity EM field, i.e., for relativistic 
photoelectron- according to which the rescattering will be negligible small in relativistic domain with the increasing 
the wave intensity p|- is not true, specially for Coulomb field, as we have obtained significant contribution even in 
the Born approximation when the impact of scattering potential is the smallest. Indeed, beyond the scope of the 
Born approximation the contribution of rescattering effect in ATI process will be more considerable (for instant, in 
the above considered GEA for scattering potential). 

In the scope of current approximation ^ ^ Za/137 the explicit ana lytic formulas for total ionization rate can be 
obtained utilizing properties of Bessel function. At the condition (3^) the argument of Bessel function X{N) >> 1 
and always X < N. Therefore the terms with A^ >> 1 and N ^ X give the main contribution in the sum ([f.]]). 
Besides, in this limit one can replace the summation over N with integration and approximate the Bessel function by 
Airy one 



^ ( - 



1/3 



Ai 



2/3 



1 - 



Af2 



(4.3) 



Turning to spherical coordinates, we carry out the integration over the Lp since there is azimutal symmetry with 
respect to the direction k (the OZ axis) and for ionization rate we have 



i6d6 



N=No 



2^2/3 {N -2Zf(k-Ilf 

N . 



it 



Ai' [y{N,9)] 



X + 



2 (A^ - 2Z) {k ■ n) 



(4.4) 



where 



2/3 



1 - 



" I fe-n 



Af2 



(4.5) 
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The y{N, 9) has a minimum as a function of A'^ and 9 and since Airy function exponentially decreases with increasing 
argument one can use Laplace method ( method of the steepest descent) in order to carry out the integration as over 
N as well as over 9. The extremum points of the function y{N, 9), i.e., the most probable values of N and 9 are 



EqUJ 



cos f „,, = 



l!{N„ 



no(iV™) 



and 



Urn = V{N„i,9,n) 



Fat 

2Fo 



2/3 



(4.6) 



(4.7) 



where Fg and Fat = Z^m?e^ are wave and atomic electric field strengths. At TV = and 9 = 9m have a peak 
for angular and energetic distribution. Let us note that the contribution of the rescattering effect to the angular 
distribution of the photoelectrons is nonessential. 

For << 1 when the wave electric field strength much exceeds the atomic one: _Fo >> Fat the main contribution 
in the integral give the areas 



59 ~ (iV„/2)-i/VVl + e'and 6N ~ 2(7V„/2)2/3 



(4.8) 



(angular and energetic widths of the peak) and for ionization rate we have an explicit formula which expresses directly 
the dependence upon the wave intensity 



34/3r2(2/3)^'^UBy V^O 



(4.9) 



For ?/m >> 1 or Fq << Fat (so called tunneling regime of ionization) we shall use the following asymptotic formula 
for Airy function 



Ai{x) 



1 



2V^ 



1/4 

X ' exp 



2x3/2 



and applying Laplace method we have 



exp 



2Fat^ 

'3 Fo 



(4.10) 



(4.11) 



Let's revert to the Born conditio n (|3.3| ) to substantia te th e condition (3.4). As is shown above we have a peak for 
angular and energetic distribution ( [l.l| ) at and Nm ( |4.6| ), and the electron mean velocity will be defined by these 
values 



rt(iv™) 



(4.12) 



■Substi tutin g ( 4.12 ) into Eq.(p.3|) we have the condition of the Born approximation in ATI process of hydrogen- like 
atoms (3.4). 

Using the explicit analytic formulas for total ionization rate we can conc lude that at = Nm and 9 = 9m have 
a peak for angular and energetic d istr ibution which are given by Eq. ( |4.6| ) with the angular and energetic widths of 
the peak 69 and 6N, respectively (E^). 



In the case of linear polarization of the wave from Eq. (3.11) we have 
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(4.13) 



where Jn{u,v) is the real generalized Bessel function (e.g., see [|18[). As is seen from the formula ( 4.13 ), in this case 
the total probability of ATI process includes all intermediate transitions of photoelectron trough the virtual vacuum 
states as well, corresponding the emission and absorption of wave photons of number — oo < n < oo (the sum over n) 
in accordance with the above mentioned behavior of wave intensity effect at linear polarization (strongly unharmonic 
oscillatory motion of photoelectron). 

Let us now consider the ATI process with the rescattering effect in nonrelativistic limit since the theoretical treat- 
ments of this problem- the main of those are Keldysh-Faisal-Reiss ansatz ( |lO]- 18|) - in general have been carried 



out for nonrelativistic photoelectron when the rescattering effect is neglected. In the pioneer result of Keldysh ( [ |16[ ) 
the rescattering of photoelectron on the potential of atomic remainder has been approximately estimated and putted 
in the form of coefficient in the ultimate formula for the ionization probability (for the wave fields much smaller than 
atomic ones). Further the same approach has been made in [ p3[ for relatively large wave fields up to the atomic ones. 
Besides, in the existing nonrelativistic theory of ATI the gauge problem for description of interaction with the wave 
field and different views concerning the role of wave intensity in the dipole approximation have arose. For discussion 
of these problems the special paper has been devoted ( |^). Moreover, in the scope of the same Keldysh-Faisal-Reiss 
ansatz the existence of stabilization effect depends on the gauge of the wave field ( |^). So that, we shall consider the 
results of the present paper in the nonrelativistic limit taking into account the photoelectron rescattering. 

From the formula ( |0 ) for the differential probability of ATI transition rate in the case of circular polarization of 
EM wave in the nonrelativistic limit we have 
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where 



(4.14) 



(4.15) 



z^ Z = Zi = e'^Al/Amuj, and Nq = {(-f^/2m - Eb) juj + z). 

The corresponding condition of the Born approximation (Eq. ( ^.4| )) in nonrelativistic limit is 



l>e> 



Z^ 
137' 



(4.16) 



The first term in the quadratic brackets of Eq. ( 4.14 ) coincides with the above-threshold ionization differential 
probability obtained in the SFA for the nonrelativistic photoelectron [|l8| without the rescattering effect. According 

to the 1^ SFA is expected to become valid when the ponderomotive potential Up ~ A^^l2m, due to EM radiation field 
larger than the ionization potential of the atom: C/p 3> Eb and consequently"!?^ /2m ^ Eb which is the condition 
o f the Born approximation. Then taking into account the scattering potential by perturbation theory we obtain 
[(4.14)] that the contribution of the photoelectron rescattering in the ATI probability (in the first order of the Born 
approximation over the Coulomb potential) is of the order of the main results of KFR ansatz. So, the neglecting of 
SB process for the photoelectron in the Coulomb field of an at omic remainder is incorrect. 

In the case of linear polarized EM wave from the formula ( 4.13 ) we have the differential probability of the ATI 
process in the nonrelativistic domain 



d^ 



Slo 



Eb 

LU 



5/2 



E 



{N- 



Eb/lo 



.1/2 



-J: 



N 



11 



{N-z-Eb/oj) 
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where 



Eb 



1/2 



COS 9 



(4.17) 



(4.18) 



in the polar coordinate system has been chosen for this case, 6 is the angle between the velocity vector of the emitted 
photoelectron and the wave polarization vector. 



APPENDIX A: DEFINITION OF THE FUNCTION Jn{U,V,A) 
A function Jn{u,v, A) may be defined by 

J„(m,w, A) = (27r)"i / decxp[i{usm{0 + A) + vsm2e -n{e + A))] 

J — TT 

or by an infinite series representation 

oo 

J„(u,«,A)= e-^''^J„-2fe(«)Jfe(t'). 



We are perform two important theorems, which can be proved from Eq. (Al): 

oo 

e"(^+'^) J„(w, V, A) = exp {i [u sin((^ + A) + ^; sin 2^]} 

n— — OO 

and 

oo 

JnTk(u,V, A)Jk{u' ,v' ,±A) = Jn{u±u',V±v',A). 



An integration by parts in Eq. (jAl|) yields to the relation 

2nJn{u,v, A) = u [J„_i(u, u, A) + J„+i(u, u, A)] 

+ 2v [e~'^^J,,_2iu, V, A) + e'^^Jn+2iu, v, A)l . 



From either Eq. (|Al|) or ( |A2[ ) follows that 



J„(m, 0, A) = J„(m), 



and 



J„(0,i;,A) 



(w), n even 
0, n odd 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 
(A6) 

(A7) 
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